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SUMMARY 

A  new  analysis  of  lunar  occultation  observations  from  1955  to  1974  utilizing 
Atomic  Time  gives  a  value  for  the  empirical  part  of  the  secular  acceleration 
of  the  Moon’s  mean  longitude  of  (  —  65"  ±  i8")/cy2.  This  differs  from  other 
determinations  which  utilized  the  Ephemeris  Time  scale,  based  on  the 
apparent  annual  motion  of  the  Sun  about  the  Earth.  These  latter,  which  give 
(  — 38"±4")/cy2,  measure  only  the  tidal  component  of  the  Moon’s  longitude 
acceleration.  The  remaining  acceleration  (  —  27"  ±  i8")/cy2,  has  as  its  most 
probable  cause  a  decrease  in  the  Universal  Gravitational  Constant  at  the  rate 
of  G/G  —  (  —  8±5)x  io-11/yr.  There  is  a  large  body  of  supporting  evidence 
for  a  rate  of  about  this  size,  including  the  disappearance  of  the  artificial  ‘  non- 
tidal  ’  component  in  the  deceleration  of  the  Earth’s  rotation,  geophysical 
evidence  for  an  observed  expansion  of  the  Earth’s  radius,  and  generally  good 
agreement  with  recent  determination  of  the  Hubble  constant,  which  is  related 
to  CjtJG.  The  observed  rate  is  also  consistent  with  the  Dirac  and  the  Hoyle- 
Narlikar  cosmological  theories,  and  to  a  lesser  degree,  with  the  Brans-Dicke 
theory.  Other  possible  interpretations  of  the  observed  excess  (negative) 
lunar  acceleration  are  also  discussed;  but  only  those  with  cosmological 
significance  seem  at  all  plausible. 

INTRODUCTION 

No  variation  in  a  fundamental  constant  of  physics  has  yet  been  demonstrated. 
However,  the  expansion  of  the  Universe  has  suggested  to  many  theoreticians  that 
the  Universal  Gravitational  Constant,  G,  may  be  decreasing,  either  as  a  cause  or 
an  effect  of  the  expansion.  Indeed,  at  least  three  currently  plausible  cosmological 
theories,  Brans-Dicke  (1961),  Hoyle-Narlikar  (1972),  and  Dirac  (1973),  demand 
that  G  be  decreasing.  In  the  strictest  classical  form  of  the  widely-accepted  general 
relativity  theory,  however,  G  must  remain  constant.  This  paper  deals  with  observa¬ 
tional  evidence  that  G  is,  in  fact,  decreasing  at  a  rate  of  about  one  part  in  io10  per 
year,  and  summarizes  some  implications  of  this  result,  along  with  other  possible 
interpretations  of  the  observational  data. 

PRINCIPLE  BEHIND  DETERMINATION 

An  immediate  consequence  of  a  decreasing  gravitational  constant  is  a  slow, 
adiabatic  expansion  of  all  orbits  in  the  solar  system  about  their  primaries.  Yinti 
(1974)  has  shown  that,  for  two-body  motion  with  G  varying  inversely  with  time 
and  averaging  over  a  complete  revolution, 

G_a_ih__iP  ,  . 

G  a  2  n  2  P’ 

where  a  dot  denotes  the  time  derivative,  a  =  semi-major  axis,  n  —  mean  motion, 

P  =  period  of  revolution.  Hoyle  (1972)  has  shown  how  such  a  change  in  the  mean 
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motion  of  the  Moon  about  the  Earth  can  be  distinguished  from  similar  changes  due 
to  other  causes,  such  as  tidal  friction.  Since  equation  (x)  holds  for  the  Earth’s 
orbit  around  the  Sun,  as  well  as  for  the  Moon’s  orbit  around  the  Earth,  it  follows 
that  an  astronomical  time  scale  derived  from  the  observed  motion  of  the  Sun  about 
the  Earth,  such  as  Ephemeris  Time,  will  suffer  a  slowing  at  the  same  proportional 
rate,  G/G.  Hence,  when  we  measure  the  acceleration  of  the  Moon’s  mean  longitude 
utilizing  Atomic  Time  (presumed  uniform),  we  obtain  the  total  acceleration  from  all 
causes.  But  when  we  measure  the  same  acceleration  utilizing  Ephemeris  Time,  any 
contribution  due  to  G  would  be  excluded,  since  it  would  be  absorbed  into  the  time 
scale.  The  difference  between  the  two  acceleration  values  would  then  indicate  the 
amount  due  exclusively  to  a  decreasing  gravitational  constant. 


OBSERVATIONAL  DATA  UTILIZING  EPHEMERIS  TIME 

There  have  been  five  determinations  of  the  acceleration  of  the  Moon’s  mean 
longitude  utilizing  Ephemeris  Time.  The  first  is  due  to  Spencer-Jones  (1939),  who 
analysed  observations  of  the  Sun,  Moon,  Mercury,  and  Venus  over  the  last  three 
centuries.  His  results  were  converted  to  their  Ephemeris  Time  equivalent  by 
Clemence  (1948),  and  the  mean  error  was  estimated  by  Morrison  (1972a),  obtaining 
a  negative  lunar  acceleration  of  ( —  22"  ±  7")/cy2  (cy  =  century),  over  and  above  the 
dynamical  acceleration  induced  by  the  slowly  decreasing  eccentricity  of  the  Earth’s 
orbit  around  the  Sun.  The  second  determination  is  due  to  Murray  (1957),  who 
analysed  ancient  observations  and  obtained  (  —  42"  ±  6")/cy2.  The  third  determina¬ 
tion  is  Newton’s  (1969),  who  analysed  ancient  observations  at  two  different  epochs, 
primarily  of  solar  eclipses,  to  obtain  independent  estimates  of  (  — 4i"-6±4"*3)/cy2 
and  (  —  42"’3  ±  6"-i)/cy2.  The  fourth  determination  was  made  by  Oesterwinter  & 
Cohen  (1972),  who  analysed  meridian  circle  observations  of  the  Sun,  Moon  and 
planets  since  1913  to  obtain  (  —  38"  +  8")/cy2.*  The  fifth  determination  is  the  recent 
result  of  Muller  &  Stephenson  (1974),  who  utilized  more  extensive  ancient  eclipse 
data  than  Newton  or  Murray,  and  used  an  improved  method  of  analysis.  Their 
result  is  (-37"- 5  ±  5"-o)/cy2. 

Of  these  five  determinations,  the  most  suspect  is  the  Spencer-Jones  value, 
which,  as  has  been  pointed  out  by  Morrison  (1972b),  is  influenced  greatly  by 
seventeenth  century  observations  of  transits  of  Mercury  across  the  Sun.  Omitting 
this  first  determination  because  of  the  high  probability  of  systematic  error,  and 
assuming  that  the  Muller  &  Stephenson  result  obsoletes  the  earlier  results  from 
ancient  observations,  the  average  of  these  last  two  determinations  gives 

«met  =  ( -  38"  ±  4")/cy2  (2) 

as  the  acceleration  of  the  Moon’s  mean  longitude,  utilizing  Ephemeris  Time.  All 
of  the  determinations,  except  that  of  Spencer-Jones,  lie  within  their  own  mean 
errors  of  this  adopted  value.  As  will  become  evident  later,  adopting  any  value 
closer  to  the  Spencer-Jones  result  would  make  the  derived  rate  of  decrease  of  G 
larger. 

*  The  statement  in  the  Oesterwinter  &  Cohen  paper  that  the  Atomic  Time  scale  was 
extrapolated  backwards  to  1913  is  incorrect  in  the  context  of  the  present  discussion — if  G 
changes,  then  the  time  scale  which  was  extrapolated  backwards  was  an  Ephemeris  Time 
scale,  adjusted  in  epoch  and  rate  to  fit  the  Atomic  Time  scale  from  1955  to  1969.  Fig.  8  of 
their  paper  should  be  re-labelled,  for  present  purposes,  as  ‘  New  Ephemeris  Time  minus 
old  Ephemeris  Time,  .  .  .’. 
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|  OBSERVATIONAL  DATA  UTILIZING  ATOMIC  TIME 

I 

I  The  analysis  of  lunar  occultation  data  utilizing  Atomic  Time  to  determine  the 
;  Moon’s  mean  longitude  acceleration  has  been  discussed  earlier  by  Van  Flandern 
;  (1970),  obtaining  (  — 52"!  i6")/cy2,  and  by  Morrison  (1973),  who  obtained 
(  — 42"  ±  6")/cy2.  The  results  of  a  new  analysis,  presented  here,  differ  from  both 
earlier  results  primarily  because  of  the  use  of  a  numerically-integrated  lunar 
ephemeris  (see  Appendix),  instead  of  the  Brown-Eckert  analytical  theory;  and 
because  of  the  inclusion  of  new  observations. 

The  occultation  observations  in  the  new  analysis  cover  the  period  1955-1974, 
during  which  Atomic  Time  was  available.  To  guard  against  an  obvious  source  of 
systematic  error,  a  solution  was  performed  only  for  photoelectric  timings  (1786 
observations),  because  insufficient  information  was  available  to  us  to  remove 
completely  the  effects  of  the  reaction  delay  of  the  observer  for  visual  timings 
before  1967.  Also,  the  effect  of  the  non-uniform  time  distribution  of  the  observa¬ 
tions  was  minimized  by  a  simple  weighting  scheme,  such  that  the  total  weight  of 
all  observations  within  each  year  was  kept  nearly  the  same.  The  complete  solution 
contained  21  unknowns,  including  selected  lunar  and  solar  elements,  lunar  size 
and  shape  parameters,  and  star  reference  system  corrections,  for  which  correlations 
were  kept  to  a  minimum  (Van  Flandern  1971).  The  total  acceleration  not  accounted 
for  by  ordinary  gravitational  theory  is  determined  to  be  (  —  65"!  io")/cy2.  The 
formal  mean  error  is  given.  This  result  was  subjected  to  significance  tests  and 
vulnerability  to  various  known  sources  of  systematic  error.  Some  of  these  are 
discussed  in  the  Appendix.  As  a  result  of  such  tests,  a  confidence  level  more  realistic 
than  the  formal  mean  error  was  determined,  giving 

«mat  =  (  —  65"+  i8")/cy2  (3) 

for  the  acceleration  utilizing  Atomic  Time. 


THE  RATE  OF  CHANGE  OF  G 

Differencing  the  lunar  acceleration  utilizing  Atomic  Time,  equation  (3),  and 
the  acceleration  utilizing  Ephemeris  Time,  equation  (2),  we  obtain  the  acceleration 
presumed  due  to  a  changing  gravitational  constant,  that  is  (  — 27"!  i9")/cy2. 
Taking  this  result,  call  it  together  with  the  lunar  mean  motion,  «m  =  x7,,*33 
x  io8/cy,  we  obtain 

-  ~~  =  ( -  o-8  ±  0-5)  x  io-8/cy  =  ( -  8  ±  5)  x  io-n/yr.  (4) 

2  n  m 

Finally,  by  equation  (1),  we  have  the  result  that  GjG  —  J(«mg/»m)-  In  what  follows, 
we  will  use  the  symbol  A  for  |(kmg/«h),  the  observed  numerical  rate  given  by 
equation  (4). 


EXAMINATION  OF  ALTERNATIVE  HYPOTHESES 

In  the  next  section,  we  will  discuss  evidence  that  a  secular  decrease  in  G  is  a 
reasonable  hypothesis,  and  is  in  fact  consistent  with  other  observational  data.  In 
this  section,  we  will  consider  the  alternative  possibilities  for  explaining  the  observed 
excess  (negative)  lunar  acceleration  when  Atomic  Time  is  utilized. 
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Examining  the  variation  of  Kepler’s  third  law, 


n  a 

2-+3 - = 

n  a 


G  M 

G+M’ 


(5) 


where  M  is  the  mass  of  the  two-body  system  being  considered,  we  see  that  there 
are  several  possible  interpretations  of  an  observed  non-zero  value  for  n/n.  The 
interpretation  preferred  here  with  G/G  ^  o  and  1VI/M  =  o  is  given  by  equation 
(x),  and  will  subsequently  be  referred  to  as  ‘  Hypothesis  G  However,  there  are 
other  interpretations  possible,  including  some  for  which  G\G  remains  zero. 

Let  us  first  consider  the  case  for  which  G/G  =  o  and  MjM  —  A.  If  all  masses 
decrease  at  the  same  rate  A,  this  situation  is  not  meaningfully  different  from 
decreasing  G.  Indeed,  all  celestial  mechanics  experiments  measure  only  the 
product  GM,  and  cannot  separate  the  two.  If  only  the  mass  of  the  Earth  (or  Earth 
plus  Moon)  were  decreasing  at  the  rate  A,  then  the  acceleration  given  by  equations 
(2)  and  (3)  should  be  equal,  which  is  clearly  not  the  case.  The  most  interesting 
possibility  for  mass  variation  is  that  only  the  mass  of  the  Sun  decreases  at  the  rate  A. 
Let  us  call  this  ‘  Hypothesis  M  (mass)  ’.  In  this  case,  the  orbits  of  all  the  major 
planets  expand,  and  their  periods  increase,  as  before;  and  therefore,  Ephemeris 
Time  will  decelerate  with  respect  to  Atomic  Time.  However,  the  orbit  of  the 
Moon  about  the  Earth  would  not  expand,  since  the  Earth’s  mass  would  not  be 
changing  significantly.  Under  this  assumption,  all  of  the  acceleration  of  the  Moon’s 
longitude  in  equation  (3)  must  be  ascribed  to  tidal  friction,  or  some  such  mechanism; 
and  the  acceleration  in  equation  (2)  is  smaller  because  of  the  slowing  of  the 
Ephemeris  Time  scale.  This  hypothesis  has  the  attractive  feature  that  it  would 
explain  the  discordance  of  the  Spencer-Jones  value  for  the  lunar  acceleration  to  the 
extent  that  it  is  actually  dependent  upon  17th  century  transits  of  Mercury,  as 
discussed  by  Morrison  (1973).  To  test  Hypothesis  M,  consider  the  possible  ways 
in  which  the  Sun’s  mass  could  be  decreasing.  From  E  —  me 2  and  the  solar  constant, 
we  can  compute  that  the  rate  of  conversion  of  matter  to  energy  in  the  Sun  causes 
a  mass  loss  of  about  7  x  io~14  solar  masses  per  year.  (A  solar  mass  is  2X  io33  g.) 
Assuming  that  the  solar  wind  is  isotropic,  the  mass  transported  away  from  the  Sun 
by  the  wind  is  at  most  about  3  x  io~13  solar  masses  per  year.  Using  recent  Sky  lab 
results  that  solar  flares  which  remove  io~17  solar  masses  from  the  Sun  per  event 
occur  with  a  frequency  of  one  per  100  hr  (Newkirk  1974,  private  communication), 
the  loss  rate  might  be  io~15  solar  masses  per  year.  Even  the  largest  solar  promi¬ 
nences  on  record  are  unlikely  to  remove  more  than  about  io-16  solar  masses. 
Therefore,  the  requirement  of  a  mass  loss  rate  of  Ax  io~10/yr  for  the  Sun  would 
seemingly  rule  out  Hypothesis  M. 

A  second  possible  interpretation,  in  which  G/G  =  o  and  M jM  =  0,  is 
dja  —  —§A.  Call  this  ‘  Hypothesis  S  (space)  ’.  This  implies  a  uniform  expansion 
of  all  space,  including  solid  bodies,  at  the  rate  of  —  §  A.  For  orbiting  bodies  only, 
the  expansion  rate  for  the  semi-major  axis  is  twice  the  radial  expansion  rate  because 
of  the  increase  in  angular  momentum.  If  we  change  the  units  of  —  §  A  for  compari¬ 
son  with  the  Hubble  constant,  we  get  52  km-1  s-1  Mpc-1.  (The  observed  Hubble 
constant  is  apparently  between  40  and  100  km-1  s_1  Mpc-1.)  In  support  of  this 
hypothesis,  Wesson  (1973)  has  summarized  observational  data  which  suggest  that 
the  Earth  has  been  expanding  at  about  the  Hubble  rate  throughout  its  existence. 
The  hypothesis  is  made  more  attractive  by  the  prediction  that  the  Earth  would 
originally  have  had  only  half  its  present  radius,  and  that  all  of  the  present  continents 
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would  then  have  formed  a  covering  over  the  entire  Earth.  Sea-floor  spreading  is 
readily  understandable  with  such  an  expanding  model  for  the  Earth.  To  test 
Hypothesis  S,  we  need  a  direct  measure  of  the  rate  of  expansion  of  some  orbit  in 
the  solar  system.  The  laser  ranging  project  for  the  Moon  will  probably  make  such 
a  measurement  before  radar  ranging  to  the  planets  does.  But  at  the  present  time, 
we  cannot  observationally  distinguish  between  Hypotheses  G  and  S.  The  principal 
reasons  for  preferring  Hypothesis  G  are  the  lack  of  a  comprehensive  cosmological 
theory  incorporating  Hypothesis  S,  and  the  continuous  increase  in  angular 
momentum  of  all  orbiting  and  rotating  bodies  which  it  implies. 

A  third  possibility  is  that  the  length  of  the  Atomic  second  is  changing  relative 
to  uniform  time.  Call  this  *  Hypothesis  T  (time)  ’.  In  essence,  we  are  postulating 
a  secular  increase  in  the  number  of  oscillations  of  the  caesium  atoms  per  Ephemeris 
second  (here,  of  constant  length).  This  would  seem  to  require  a  variation  in  some 
fundamental  constant  of  physics  other  than  G.  Moreover,  it  is  probably  immaterial 
whether  we  consider  the  Atomic  second  as  constant  and  other  physical  constants 
as  varying,  or  vice  versa. 

Dicke  (private  communication)  has  pointed  out  that  variable  tidal  friction 
is  another  possible  explanation — ‘  Hypothesis  F  ’.  This  is  quite  true.  However, 
variation  by  a  factor  of  50  per  cent  within  a  few  decades  is  not  likely.  This  can  be 
stated  with  even  greater  confidence  in  the  light  of  the  new  Muller  &  Stephenson 
results,  which  show  that  the  tidal  deceleration  has  remained  essentially  constant 
for  over  3000  yr.  Any  extremely  short  term  variations  should  average  out  over  the 
Moon’s  nodal  period  of  18-6  yr. 

There  are  an  unlimited  number  of  possible  ways  in  which  Kepler’s  law  can  be 
violated,  due  to  the  action  of  small  forces  which  have  not  been  considered  in  the 
equations  of  motion.  In  this  category  we  might  consider  solar  radiation  pressure, 
solar  wind,  magnetic  forces,  meteoritic  impacts,  passing  comets,  the  tidal  action 
of  the  Sun  on  the  Earth,  and  undiscovered  masses  in  the  solar  system.  In  each  of 
the  cases  just  mentioned,  estimates  of  the  reasonable  upper  limit  for  the  size  of 
these  effects  are  negligible.  None  the  less,  the  main  reason  for  preferring  Hypothesis 
G  to  the  ‘  unmodelled  force  ’  hypothesis  is  the  supporting  evidence  which  follows. 


DISCUSSION  OF  SUPPORTING  EVIDENCE 

Hoyle  (1972)  has  reviewed  much  of  the  evidence  in  support  of  Hypothesis  G, 
that  the  gravitational  constant  is  decreasing  at  about  a  part  in  io10  per  year.  Perhaps 
the  most  persuasive  implication  is  that  of  conservation  of  angular  momentum  in 
the  Earth-Moon  system.  If  we  take  Muller  &  Stephenson’s  (1974)  results,  for 
example,  then  io9(co/to)  for  the  Earth’s  rotation  is  about  (  — 29±3)/cy,  when 
measured  utilizing  Ephemeris  Time.  (Newton’s  (1973)  result  was  —2513.)  This 
breaks  down  into  a  frictional  part  of  (  —  44  +  5)/cy  from  the  tidal  action  of  the  Sun 
and  Moon,  and  a  non-frictional  part  of  (+i5±6)/cy  of  unknown  origin.  (The 
frictional  part  comes  from  equation  (2),  by  applying  Newton’s  relation 

io9(cb/w)  =  1-165  ^met, 

as  required  to  conserve  angular  momentum.)  When  Atomic  Time  is  utilized 
instead,  observed  values  of  cb/co  utilizing  Ephemeris  Time  must  be  corrected  by 
2 A,  where  A  is  given  by  equation  (4).  (The  slowing  of  the  Earth’s  rotation  due  to 
the  expansion  which  accompanies  a  decreasing  G  is  a  much  smaller  effect.)  Hence, 
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Muller  &  Stephenson’s  results  for  io9(cb/a>)  would  change  from  (  —  29  +  3)/cy  to 
(  — 45  +  io)/cy,  in  excellent  agreement  with  the  frictional  part  of  (  —  44+  5)/cy.  No 
unexplained  ‘  non-frictional  ’  change  in  the  Earth’s  rotation  remains. 

Other  arguments  cited  by  Hoyle  in  support  of  Hypothesis  G  include  the 
thermal  history  of  the  Earth,  which  would  have  been  much  hotter  in  its  early 
history  than  at  present  (consistent  with  the  evolutionary  trend  of  maximum 
survival  temperatures  for  the  most  primitive  life  forms  on  Earth);  the  calculation 
that  the  initial  helium  content  of  the  Sun  would  have  to  be  about  15  per  cent 
rather  than  about  25  per  cent  to  reach  its  present-day  luminosity  (consistent  with 
solar  wind  measures  that  the  outer  solar  layers,  which  are  believed  to  have  remained 
essentially  unchanged  over  the  lifetime  of  the  solar  system,  contain  15  per  cent 
helium);  and  the  prediction  of  a  blue,  starlike  class  of  quasi-stellar  galaxies  (such 
as  recently  found  by  Sandage)  due  to  an  enhancement  of  the  rate  of  stellar  evolution 
as  a  function  of  mass. 

Wesson  (1973)  cites  further  geophysical  evidence,  including  an  illuminating 
discussion  of  continental  drift  and  sea-floor  spreading.  If  G  is  decreasing,  the  Earth 
must  expand  somewhat  due  to  the  decreased  weight  of  the  surface  layers,  but  at 
perhaps  only  10-20  per  cent  of  the  rate  in  Hypothesis  S.  Wesson  points  out  that 
the  expanding  globe  model  makes  it  easier  to  understand  how  a  mid-oceanic  ridge 
system  can  almost  surround  a  continent  like  Antarctica. 

Morganstern  (1972)  has  examined  a  pressure-filled  curved-space  cosmological 
situation,  and  finds  that  positively-curved  spaces  with  pressure  predict  a  GjG  rate 
in  agreement  with  equation  (4). 

Newton  (1968)  also  came  to  the  conclusion  that  the  Earth’s  polar  moment  of 
inertia  is  increasing  and  that  G  is  decreasing  at  rates  quite  consistent  with  equation 
(4),  and  with  the  rate  of  upward  transport  of  material  at  the  mid-ocean  ridges. 

It  should  also  be  mentioned  that  if  Dicke  (1966)  had  used  the  modern  data  for 
«met,  instead  of  the  Spencer-Jones  value,  in  his  study  of  the  Earth’s  rotational 
history,  his  estimated  G/G  values  would  have  been  in  good  agreement  with  the 
present  result,  instead  of  the  smaller  value  suggested  by  the  Scalar-tensor  theory 
(Brans  &  Dicke  1961). 

Cosmological  theories  consistent  with  Hypothesis  G  are  those  of  Hoyle- 
Narlikar  (1972),  Dirac  (1973),  and  Brans-Dicke  (1961).  These  all  have  the  important 
features  that  they  are  nearly  consistent  with  general  relativity  in  so  far  as  the  solar 
light  bending  and  Mercury  perihelion  advance  tests  are  concerned;  and  they  are 
characterized  by  a  properly  integrated  form  of  Mach’s  principle.  In  the  Hoyle- 
Narlikar  theory,  equation  (4)  gives  twice  the  Hubble  expansion  rate.  Therefore, 
the  Hubble  constant  implied  by  equation  (4)  would  be  (39  ±  24)  km-1  s-1  Mpc-1. 
For  the  other  two  cosmologies,  the  Hubble  constant  would  equal  the  solar  system 
expansion  rate  of  (78149)  km-1  s-1  Mpc-1,  according  to  equation  (1).  These  are 
in  reasonable  agreement  with  some  of  the  most  recent  determinations  of  that 
parameter.  Sandage  (1972),  for  example,  has  obtained  (55  +  7)  km-1  s_1  Mpc-1. 
The  Hubble  age  of  the  Universe,  from  equation  (4),  would  be  (12!  41)  x  IC)9  Yr- 


FUTURE  OBSERVATIONAL  DATA 

In  the  solar  system,  many  of  the  observational  consequences  of  a  decreasing  G 
increase  with  time  squared.  In  addition,  correlations  between  accelerations  and 
other  parameters  decrease  rapidly  with  time.  In  the  case  of  lunar  occultation  data, 
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|  the  time  base  begins  in  1955,  when  Atomic  Time  came  into  use.  About  two  more 
I  years  of  data  will  result  in  halving  the  mean  error  of  the  present  determination. 

Independent  determinations  are  now  nearly  possible  from  planetary  radar 
;  ranging  results.  A  determination  by  Shapiro  et  al.  (1971)  set  an  upper  limit  of 
:  GjG  4  X  io-10/yr.  A  more  recent  estimate  (Shapiro  1974)  from  Mercury  ranging 
had  decreased  this  estimate  to  (  +  4±8)x  io-11/yr,  but  with  the  separation  of  the 
acceleration  from  the  topography  of  Mercury  not  yet  complete.  The  lunar  laser 
ranging  program  is  also  very  close  to  obtaining  a  lunar  acceleration  measure,  and 
should  have  preliminary  results  by  the  end  of  1974.  Moreover,  some  laboratory 
experiments  are  now  contemplated  for  the  next  few  years  which  can  reach  better 
than  1  part  in  io137yr. 

On  the  theoretical  front,  the  consequences  of  G  varying  must  be  reviewed 
carefully  in  the  fields  of  geophysics,  astrophysics  and  cosmology  for  possible 
contradiction  with  observational  data,  as  in  the  discussion  of  the  Hoyle-Narlikar 
theory  by  Barnothy  &  Tinsley  (1973).  Since  the  possibility  of  G  decreasing  has 
now  become  of  significant  probability,  it  will  be  important  to  make  the  effort  to 
incorporate  varying  G  into  stellar  evolution  calculations,  for  example,  as  was  done 
by  Pochoda  &  Schwarzschild  (1964).  If  the  G  Hypothesis  is  contradicted,  then 
alternate  hypotheses,  especially  Hypothesis  S,  should  be  tested. 

CONCLUSION 

If  Atomic  Time  is  a  uniform  time  scale,  then  the  observed  discordance  between 
the  acceleration  of  the  Moon’s  mean  longitude  utilizing  Ephemeris  Time  with  the 
same  acceleration  utilizing  Atomic  Time  necessarily  implies  a  deceleration  in  the 
Ephemeris  Time  scale;  or  equivalently,  there  must  be  an  unmodelled  acceleration 
of  the  Sun’s  apparent  annual  motion  about  the  Earth.  If  Ephemeris  Time,  Te, 
and  Atomic  Time,  Ta,  are  adjusted  so  as  to  be  equal  in  epoch  and  rate  at  some 
moment  To  (with  time,  T,  in  centuries),  then 

7a  =  A(T—  T0)2,  (6) 

where  A  is  given  by  equation  (4).  If  we  wish  to  express  the  difference 
(Te—  Ta)  =  ATa  in  seconds,  we  would  have 

ATa  =  —  25S-25(T— T0)2.  (7) 

The  corresponding  acceleration  in  the  Sun’s  mean  longitude  utilizing  Atomic 
Time  is  given  by 

wsAT  =  —  2"-07/cy2,  (8) 

«sET  being  zero  by  definition.  The  evidence  presented  here  suggests  that  the  most 
probable  cause  of  the  anomalous  accelerations  is  a  secular  decrease  in  the  Universal 
Gravitational  Constant,  G,  at  the  rate  of  (  — 8±5)x  io_11/yr.  This  interpretation 
of  the  cause  is  not  compelling;  but  it  seems  to  be  strongly  indicated,  since  the 
anomalous  accelerations  almost  certainly  have  a  cosmological  cause,  and  since 
there  is  other  supporting  evidence. 
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APPENDIX 

This  analysis  utilizes  a  numerically-integrated  lunar  ephemeris  in  the  com¬ 
parison  with  occultation  observations.  The  integrated  ephemeris  originally  adopted 
was  that  of  Garthwaite,  Holdridge  &  Mulholland  (1970).  This  will  be  referred  to 
subsequently  as  LE16.  The  result  for  G/G  originally  obtained,  using  LE16,  was 
(  —  i2±3)x  io_11/yr.  It  was  pointed  out  by  a  number  of  people  that  the  validity 
of  LE16  was  of  crucial  importance  to  the  principal  result  of  this  paper.  Furthermore, 
there  was  evidence  in  the  literature  (see,  e.g.  Oesterwinter  &  Cohen  1972;  Bender 
et  al.  1973)  that  there  were  discrepancies  between  LE16  and  the  Oesterwinter  & 
Cohen  lunar  integration,  which  we  will  call  NWL.  Therefore,  a  comparison  of  the 
two  integrations  was  performed. 

The  differences  NWL-LE16  were  formed  during  the  23-yr  period  (1950-73) 
common  to  the  two  integrations.  The  effects  of  completely  different  starting  con¬ 
ditions  or  elements  for  both  the  Earth  and  the  Moon  had  to  be  removed,  along 
with  the  effects  of  differing  adopted  masses  and  other  constants,  utilizing  analytical 
partial  derivatives.  The  remaining  signal  was  still  highly  systematic.  A  power 
spectrum  analysis  was  then  performed  to  identify  the  precise  periodicities,  which 
in  turn  led  to  an  understanding  of  their  causes. 

It  was  found  that  LE16  had  used  a  force  function  for  the  perturbations  on  the 
Moon  due  to  the  oblateness  of  the  Earth  which  omitted  the  effects  of  the  finite  mass 
of  the  Moon.  This  had  been  discovered  earlier  by  Jim  Williams  at  the  Jet  Propulsion 
Laboratory  in  connection  with  shorter  numerical  integrations  used  for  the  laser 
ranging  project.  Moreover,  it  was  found  that  the  NWL  integration  had  omitted  the 
effect  of  the  precession  of  the  Earth’s  figure,  when  calculating  the  same  perturba¬ 
tions.  After  correction  for  these  two  effects,  the  remaining  differences  between  the 
two  integrations  were  all  numerically  less  than  o"-o2. 

Since  a  problem  had  arisen  with  each  of  the  integrations,  a  third  comparison 
integration  was  considered  to  be  of  significant  value.  The  only  integration  resulting 
from  the  laser  ranging  project  which  had  been  extended  over  a  sufficiently  long 
time  period  was  LE25,  which  was  kindly  supplied  by  Jim  Williams.  The  differences 
LE25-LE16  were  formed  over  the  interval  1950-75,  and  analysed  as  before.  The 
results  differed  from  the  other  two  integrations  in  the  crucial  acceleration  term. 
Consultations  with  Jim  Williams  again  led  to  the  cause — in  this  case,  a  libration 
theory  truncated  too  soon  in  LE25,  so  that  the  mean  orientation  of  the  lunar  bulge 
toward  the  Earth  was  not  exactly  in  the  mean  Earth  direction.  This  problem  was 
not  yet  of  any  significant  importance  to  the  laser  ranging  project.  After  applying 
this  correction,  the  three  integrations  mutually  agreed  to  within  o"-02. 

The  three  force  function  problems,  discovered  in  this  analysis  (one  in  each 
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integration)  resulted  in  discordances  which  looked  like  sin  £2,  cos  £2,  and  Tz, 
respectively,  where  £2  is  the  mean  longitude  of  the  Moon’s  ascending  node  on  the 
ecliptic  (period  18-6  yr).  Because  of  their  long  period,  both  £2  terms  are  correlated 
with  the  T2  term,  which  measures  the  acceleration.  Hence,  all  three  problems 
would  have  influenced  the  determination  of  GjG. 

Although  the  full  effect  of  each  of  these  known  problems  has  been  removed, 
and  the  three  integrations  now  agree  with  each  other,  it  is  clear  in  hindsight  that 
the  original  mean  error  estimate  ( ±  3  x  io-11/yr)  did  not  adequately  allow  for  such 
possible  sources  of  systematic  error.  A  mean  error  of  ( ±  5  x  io-11/yr)  would  have 
covered  each  of  the  errors  in  O/G  resulting  from  force  function  problems.  More¬ 
over,  systematic  errors  in  star  positions  also  have  the  potential  to  introduce  sin  £2 
terms,  because  the  Moon  occults  differing  declination  zones  of  stars  with  an  £2 
period.  Therefore,  based  upon  the  above  experience,  and  upon  detailed  studies  of 
possible  systematic  error  sources  in  the  observations,  the  author  believes  that 
( ±  5  x  io-11/yr)  is  a  more  realistic  mean  error  than  the  formal  solution  mean  error. 
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